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An approach to fast fits of the unintegrated gluon density has been developed and used 
to determine the unintegrated gluon density by fits to deep inelastic scatting di-jet 
data from HERA. The fitting method is based on the determination of the parameter 
dependence by help of interpolating between grid points in the parameter-observable 



The substructure of the proton is parameterized by parton density functions (PDFs). In 
'■ perturbativc QCD the PDFs are given by solutions of integral equations, for which the 

, initial input distributions have to be determined by fits. It turns out that, for exclusive 

final states, it is not statistically efficient to tune Monte Carlo event generators (MC) by 
sequential calls of the generator together with a minimisation program. Motivated by [T] , 
■ we use an alternative fitting method, which is based on producing a grid in parameter- 

r*"** , observable space. This allows the parameter dependence to be determined by polynomial 

interpolation before the fit is performed, which significantly reduce the time for performing 
the fit. 

Here we determine the parameters in the starting distribution of the unintegrated gluon 
density function (uGDF) by using the CASCADE Monte Carlo event generator [J]. The fit 
is performed to low Q 2 di-jet data from the HI experiment [3]. 



2 The Unintegrated Gluon Density 

The starting distribution of the unintegrated gluon density is parameterized as 

A (x,k t ) = Nx~ B (l - x) c cxp(k t -nf/a 2 

where x is the longitudinal momentum fraction of the proton carried by the gluon and k t its 
transverse momentum. Here we use the fitting method to determine the normalisation, N, 
the small x behaviour, B, and the shift, fi, of the Gaussian for the transverse momentum of 
the non-perturbative gluon. The parameters C and cr, is kept fixed at C = 4 and a = 2. 

The starting distribution is evolved to higher scales by gluon emissions according to the 
CCFM evolution equation which impose angular ordering of the emitted gluons. 

3 The Fitting Method 

In the first step of the fitting procedure we build up a grid of MC predictions in the parameter 
space (pi,p2, ■ ■ ■ ,Pn) for each of the observables X. 
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Secondly, we use the grid to describe the parameter-observable space analytically, by 
using a polynomial of the form 

n n 71— 1 ra 

X(pi,P2, -,Pn) 

i— 1 i— 1 i— 1 j— 2+z 

It is here possible to minimise \ 2 with singular value decomposition (SVD)[J], since the 
coefficients which are to be determined (Aq, Bi,. . . ) are mapping the MC grid on the ob- 
servables in an over determined system of linear equations. This is done separately for each 
of the MC predicted data points. 

In order to account for correlations between parameters the form of the polynomial has 
to be of at least second order. In the presented fit we use a third order polynomial, which 
gives a significantly better description of the parameter space. 

Having determined the polynomial describing the parameter space, we can fit the pa- 
rameters pi,p2, ■ ■ ■ ■ This is done by applying a x 2 minimisation to 

,poly — X k ^ data ) 2 
(3Xk,poly ~ SXk^data) 2 

where the sum runs over all data points, k. X^^data is the measured data point, with the 
corresponding experimental error 5 X^, data, and Xk tPO i y the polynomial prediction, with the 
error 8Xk tPO iy calculated from the individual errors of the fitted coefficients by using the 
covariance matrix. To perform this last step we use MINUIT since the dependence on 
parameters p\,p2, ... is non-linear. 

The method turns out to be very time efficient in particular since the MC grid points 
are generated simultaneously and instead of fitting MC to data we perform a fit of the 
polynomials, which is much faster due to the fact that the event generation already is 
performed. The method also allows for very fast refitting if one, for example, wants to study 
the exclusion of some experimental data points. 



4 The Experimental Data 

The measurement used for the fit was made by HI [3] at ^fs = 318 GcV within the kinematic 
range 5 < Q 2 < 100 GcV 2 , 10~ 4 < x Bj < 10~ 2 and 0.1 < y < 0.7. Jets are defined in the 
7P rest frame with the inclusive fc t -algorithm, and required to fulfill E^ CAI > 5 GeV and 
— 1 < rj LAB < 2.5, where E^ CM is the transverse momentum in the hadronic center of mass 
frame and r\ LAB is the pseudorapidity in the laboratory frame. Events with at least two 
jets fulfilling these requirements were analysed. The total di-jet cross-section was measured 
as a function of A which gives an additional restriction on the hardest di-jet according to 
fiHCM > 5 + A GcV. The cross-section were made double differential by binning in x Bj and 
Q 2 . 



5 Results 

The result of the fit to the HI di-jet data is shown in Fig. [1] together with the data as well as 
three already existing uGDFs. The parameter values determined from the fit are N = 0.28, 
B = 0.25 and ji = 3.0. Scanning \ 2 as a function of these parameters, as shown in Fig. [2j 
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Figure 1: Di-jet data as a function of the difference between the transverse momentum 
requirement on the di-jets, A, in bins of XBj and Q 2 compared to predictions from the 
CASCADE Monte Carlo event generator using the newly fitted PDF (full line) and 3 old 
PDFs; set AO (dashed), set C (dotted/dark), J2003 set 2 (dotted/ light) . 
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confirms that this is a minimum for N and B, and \x has flattened out at 3. Varying this 
parameter around 3 thus not change the MC predictions. 

The x 2 /ndf of the new fit is 2.1, which should be compared to x 2 /ndf =3.5 for the best 
performing old uGDF, which in this case is set AO. A uGDF that was determined from 
fits to the proton structure function [2] . There are two significant differences between these 
two uGDFs. The first one is the small x behaviour given by B, which is for set AO, but 
suggested to be 0.25 by the fit to di-jet data. The second difference is that, while set AO 
uses a non shifted Gaussian for the gluon k t , i.e. fi = 0, the di-jet data suggests a large shift 
resulting in a decreasing starting distribution towards low k t . The differences are illustrated 
in Fig. [3] where the uGDF is draw as a function of x and k 2 . 




Figure 2: x 2 profiles as a function of the fitted parameters. 




Figure 3: The newly fitted PDF (dashed line) compared to the old PDF set AO (full line), 
drawn as a function of x (for kf =10 GeV 2 ) and k 2 (for x=0.001). 
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